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ABSTRACT 

We have investigated the Cosmic Infrared Background (CIB) anisotropies in the frame- 
work of the physical evolutionary model for proto-spheroidal galaxies by Granato et al. 
(2004). After having re-calibrated the cumulative flux function dS/dz at A ^ 850 /im 
using the available determinations of the shot noise amplitude (the original model 
already correctly reproduces it at shorter wavelengths) the CIB power spectra at 
wavelengths from 250 /xm to 2 mm measured by Planck, Herschel, SPT and ACT ex- 
periments have been fitted using the halo model with only 2 free parameters, the 
minimum halo mass and the power-law index of the mean occupation function of 
satellite galaxies. The best-fit minimum halo mass is log(Mi„in/Af0) — 12.24 ± 0.06, 
higher than, but consistent within the errors, with the estimate by Amblard et al. 
(2011) and close to the estimate by Planck Collaboration (2011). The redshift evolu- 
tion of the volume emissivity of galaxies yielded by the model is found to be consistent 
with that inferred from the data. The derived effective halo mass, Mcs — 5 x 10^^ Mq, 
of z ~ 2 sub-millimeter galaxies is close to that estimated for the most efhcient star- 
formers at the same redshift. The effective bias factor and the comoving clustering 
radius at z ~ 2 yielded by the model are substantially lower than those found for a 
model whereby the star formation is fueled by steady gas accretion, but substantially 
higher than those found for a merging-driven galaxy evolution with a top-heavy initial 
mass function. 

Key words: submillimetre: galaxies - galaxies: statistics - galaxies: haloes - galaxies: 
high redshift. 



1 INTRODUCTION 

The Herschel surveys have allowed clustering studies (Mad- 
dox et al. 2010; Cooray et al. 2010) of sub-millimeter galax- 
ies with a statistics at least one order of magnitude bet- 
ter than previously possible (Blain et al. 2004; Scott et al. 
2006). These studies have been complemented by determi- 
nations of the angular power spectrum of the Cosmic In- 
frared Background (CIB) anisotropies on BLAST (Viero et 
al. 2009), Planck (Planck Collaboration 2011), and Herschel 
(Amblard et al. 2011) maps. Due to the unique power of sub- 
millimeter surveys in piercing the distant universe, thanks 
to the strongly negative K-correction, the clustering prop- 
erties contain signatures of the large scale structure at high 
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redshifts and can allow us to discriminate between differ- 
ent formation mechanisms that have been proposed for sub- 
millimeter galaxies. For example, merger driven galaxy evo- 
lution models, that follow the evolution of both the disk and 
the spheroidal components of galaxies, predict much lower 
clustering strengths for sub-mm galaxies (e.g. Almeida et al. 
2011; Kim et al. 2011) than models whereby the star forma- 
tion is fueled by steady accretion of large amounts of cold 
gas (e.g. Dave et al. 2010). 

In this paper, building on the work by Negrello et al. 
(2007), we investigate the constraints set by mm and sub- 
mm clustering data on the physical model worked out by 
Granato et al. (2001, 2004) and further elaborated by Lapi 
et al. (2006) and Mao et al. (2007). 

A specific prediction of the model is that high-2 massive 
proto-spheroidal galaxies dominate the sub-mm counts over 
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a limited flux density range (cf. Lapi et al. 2011). At 250 
the Euclidean normalized differential counts of these objects 
peak at « 30 mjy; above ~ 60 mjy and below ~ 10 mjy 
the counts are dominated by z 1.5 quiescent and star- 
bursting late-type galaxies, less massive and less clustered 
than the high-z proto-spheroidal galaxies. The flux density 
range where proto-spheroidal galaxies dominate broadens 
and the peak shifts to brighter flux densities with increasing 
(sub-)mm wavelength. Therefore, in this scenario, the ex- 
pected clustering strengths depend on the flux density range 
that is being probed and on wavelength. 

Several other analyses of data on the angular correla- 
tion function of (sub-)mm sources and of the power spec- 
trum of the CIB anisotropies have been carried out. They 
however use phenomenological parameterized models for the 
evolution of extragalactic sources (Hall et al. 2010; Planck 
CoUaboration 2011; Millea et al. 2011; Penin et al. 2011) or 
even of the clustering power (Addison et al. 2011). Also data 
at different wavelengths are usually fitted separately (Planck 
Collaboration 2011; Amblard et al. 2011). On the contrary, 
the present analysis relies on a physical model for the evo- 
lution of proto-spheroidal galaxies (although the treatment 
of spiral and starburst galaxies is phenomenological) and 
aims at accounting simultaneously for clustering data over 
a broad range of wavelengths, from 250 /im to a few mm. 

It should be noted, however, that the physical model 
is exploited only to compute the cumulative flux function 
that weights the redshift-dependent spatial power spectrum 
in the Limber approximation for the angular power spec- 
trum. The Halo Occupation Distribution (HOD), which is 
a statistical description of how dark matter halos are pop- 
ulated with galaxies, is dealt with in a simplified manner, 
without including the relationship between luminosity and 
halo mass. This is the standard practice, justified by the 
complexity of a thorough treatment that does not appear to 
be required by existing data. In the Granato et al. (2004) 
model the star-formation rate is related to the halo mass, to 
the virialization redshift and to the age of the galaxy. Includ- 
ing these additional ingredients in the analysis is impractical 
at the present stage. A pioneering model that explicitly in- 
cludes a relationship between infrared luminosity and halo 
mass has been presented by Shang et al. (2011). 

The plan of the paper is the following. In §[2] we present 
a short overview of the evolutionary model for the relevant 
galaxy populations. In §[3] we describe the halo model for- 
malism used to compute the contributions to the power spec- 
trum of Cosmic Infrared Background (CIB) anisotropies and 
to the angular correlation function of detected galaxies (§|4]). 
Our main results are presented in §[5] and our main conclu- 
sions are summarized in §[§1 

We adopt a standard flat ACDM cosmology with h — 
//o/100kms~^ Mpc~^ — 0.70 and a local matter density 
n^o = 0.27. 



2 OVERVIEW OF THE MODEL 

The sub-millimeter extragalactic sources are a mixed bag of 
various populations of dusty galaxies and of flat-spectrum 
radio sources (see, e.g., Lapi et al. 2011). 

Our model interprets powerful high-2 sub-mm galaxies 
as massive proto-spheroidal galaxies in the process of form- 



ing most of their stellar mass (see also Blain et al. 2004; 
Narayanan et al. 2010; Dave et al. 2010). It hinges upon high 
resolution numerical simulations showing that dark matter 
halos form in two stages (Zhao et al. 2003; Wang et al. 2011; 
Lapi & Cavaliere 2011). An early fast collapse of the halo 
bulk, including a few major merger events, reshuffles the 
gravitational potential and causes the dark matter and the 
stellar component to undergo (incomplete) dynamical re- 
laxation. A slow growth of the halo outskirts in the form of 
many minor mergers and diffuse accretion follows; this sec- 
ond stage has little effect on the inner potential well where 
the visible galaxy resides. 

The star formation is triggered by the fast col- 
lapse/merger phase of the halo and is controlled by self- 
regulated baryonic processes. It is driven by the rapid cool- 
ing of the gas within a region of ~ 70(Mh/10^^ M0)^/^[(1-|- 
z)/3]~^ kpc, where Alh is the halo mass, is regulated by the 
energy feedback from supernovae (SNe) and Active Galac- 
tic Nuclei (AGNs), is very soon obscured by dust and is 
stopped by quasar feedback. The AGN feedback is relevant 
especially in the most massive galaxies and is responsible 
for their shorter duration (5 — 7 x 10** yr) of the active star- 
forming phase. In less massive proto-spheroidal galaxies the 
star formation rate is mostly regulated by SN feedback and 
continues for a few Gyr. 

Since spheroidal galaxies are observed to be in passive 
evolution at z < 1 — 1.5 (e.g., Renzini 2006), they are visi- 
ble at sub-mm wavelength only at high redshifts. Lapi et al. 
(2011) have shown that the Granato et al. (2004) model, as 
further elaborated by Lapi et al. (2006), provides a reason- 
ably good fit to the observed counts from 250 fim to ~ 1 mm 
as well as to the luminosity functions in the range z — 1 — 4 
and to the redshift distributions ai z > 1 estimated from 
Herschel-ATLAS (Eales et al. 2010) data. 

The fit was obtained using of a single SED (that of the 
well studied z = 2.3 strongly lensed galaxy SMM J2135- 
0102, "The Cosmic Eyelash"; Ivison et al. 2010, Swinbank 
et al. 2010) for the whole population of proto-spheroidal 
galaxies. This is obviously an oversimpliflcation and indeed 
the Lapi et al. (2011) counts are somewhat high at mm 
wavelengths, especially at relatively bright flux densities. As 
a consequence, the model overestimates the Poisson (shot- 
noise) contribution to the power spectrum of intensity fluc- 
tuations since such contribution is directly related to the 
source counts [see eq. (|2ip ] . Consistency with the shot-noise 
levels estimated by Planck Collaboration (2011) at 353 GHz 
and measured by Hall et al. (2010), Dunkley et al. (2011) at 
220 and 150 GHz is recovered scaling down the cumulative 
flux function dS/dz [see eq. (|19|l ] of proto-spheroidal galax- 
ies by a factor of 0.81, 0.71, and 0.55 at 353, 217, and 150 
GHz (850 pim, 1.38mm, 2mm), respectively. No correction 
was applied at higher frequencies. In practice, we use the 
determination of the shot noise amplitude to recalibrate the 
function dS/dz to be used to compute the clustering power 
spectrum, which is measured independently. This correction 
mimics the effect of adopting a SED decreasing with increeis- 
ing wavelength beyond the peak a bit more steeply than the 
one adopted by Lapi et al. (2011). 

As suggested in the latter paper, the overestimate of 
mm-wave counts may be cured if higher-2 galaxies, that 
yield larger and larger contributions to the bright counts at 
increasing mm wavelengths, have SEDs slightly hotter than 
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SMM J2135-0102 and closer to that of G15.141 (Cox et al. 
2011; see Fig. 2 of Lapi et al. 2011). We have checked that 
indeed a good fit of the counts at all the frequencies consid- 
ered here is obtained using the SMM J2135-0102 SED for 
galaxies at z < 3.5 and the SED of G15.141 at higher z. 
However the match of the frequency spectrum of the shot- 
noise amplitude also improves but not enough to reach con- 
sistency with observational estimates at the longer wave- 
lengths. Since the shot noise amplitude can be computed 
directly from the counts, this suggests that there may be 
some small, but non-negligible, offsets between the calibra- 
tion of point source flux densities and that of the diffuse 
background. This is not surprising since, in addition to the 
possibility of an imperfect photometric calibration, at mm 
wavelengths the recovery of the contribution of dusty galax- 
ies to the power spectrum requires a delicate subtraction of 
the other components (Cosmic Microwave Background, cir- 
rus emission, fluctuations due to radio sources). A rescaling 
to match the shot noise spectrum seems to be the only prac- 
tical way for correcting for these offsets. Since the modifica- 
tion is only significant at ^ 850 /im (in the observer frame), 
i.e. well beyond the peak for most sources, the impact on the 
bolometric luminosity, which is related to the halo mass, is 
minor. For galaxies at redshifts up to 2: = 3.5, accounting 
for essentially all the signal, the bolometric luminosity varies 
by 2%. For comparison, the coefficient of the relationship 
between the star formation rate (SFR; given by the model) 
and the bolometric luminosity has an uncertainty of ~ 30% 
(Kennicutt 1998). 

The Granato et al. (2004) model is meant to take 
into account the star formation occurring within galactic 
dark-matter halos virialized at z^ir ^ 1.5 and bigger than 

Mvir — W^^'^Mq, which are, crudely, associated to massive 
spheroidal galaxies. We envisage disk (and irregular) galax- 
ies as associated primarily to halos virializing at z^h 1.5, 

which have incorporated, through merging processes, a large 
fraction of halos less massive than 10^^'^ Mq virializing at 
earlier times, which may become the bulges of late type 
galaxies. The model, however, does not follow the formation 
and evolution of disk and bulge components of galaxies. For 
spiral and starburst galaxies we adopt the phenomenologi- 
cal model described by Negrello et al. (2007). On the other 
hand, as shown in the following, these galaxies are essen- 
tially non influential for the purposes of the present paper 
in the considered frequency range: proto-spheroids dominate 
the contributions both to the power spectrum of fluctuations 
and to the angular correlation function of detected sources. 

Because of the strong dilution due to their very broad 
luminosity function, the contribution of radio sources to 
the clustering power spectrum can be safely neglected in 
the wavelength range considered here. Their contribution to 
Poisson fluctuations was computed using the De Zotti et al. 
(2005) model. 



3 HALO MODEL FORMALISM 

To compare the clustering properties expected from our 
model with observational data we adopt the halo model for- 
malism (Cooray & Sheth 2002). The power spectrum of the 
galaxy distribution is parameterized as the sum of the 1- 



halo term, that dominates on small scales and depends on 
the distribution of galaxies within the same halo, and the 
2-halo term, that dominates on large scales and is related to 
correlations among different halos: 
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where dn/dM is the halo mass function (Sheth & Tormen 
1999) and the linear matter power spectrum, Piin(fc,2), has 
been computed using the CAMB cod^ (Lewis, Challinor & 
Lasenby 2000). Here, Ugai(fc,M) denotes the Fourier trans- 
form of the mass density profile of the galaxy distribution 
within the dark matter halo, that we assume to be approx- 
imately the same as that of the dark matter, i.e. we take 

ttgal(fc, M) ~ Udm(fc, M). 

(A'gai) is the mean number of galaxies in a halo of 
mass M, subdivided in "central" and "satellite" galaxies 
((A^'gai) = (A^'con) + (A^sat)), whilc figai is the mean number 
density of galaxies: 

We model the HOD using a central-satellite formalism (see, 
e.g., Zheng et al. 2005); this assumes that the first galaxy 
to be hosted by a halo lies at its center, while any remain- 
ing galaxies are classified as satellites and are distributed 
in proportion to the halo mass profile. Following Tinker & 
Wetzel (2010), the mean occupation functions of central and 
satellite galaxies are parameterized as: 



(Ace„) = 
(A,at) = 



1 + erf 
1 + erf 



logio(Af/Mmin) 

a(logio M) 

logio(Af/2Mmin 

a(logio M) 



f— V 



(5) 
(6) 



where Mmin, cvsat, Maat, and a{\og^Q M) are free parameters 
assumed to be redshift independent. In this formalism halos 
below Afniin do not contain galaxies while halos above this 
threshold contain a central galaxy plus a number of satellite 
galaxies with a power-law mass function with slope asat. 

The mean mass density profile of halos of mass M is 
(Navarro, Frenk, & White 1996): 



p(r) 



(r/r-,)(l + r/r,)2' 



M = A-Kpsrt 



log(l + c) - 



l-fc 



(7) 



(8) 



with c = rvir/^a. The normalized Fourier transform of this 
profile is: 

i*dm(fc, M) = {sin(fcr,) [^^([l + c]kr,) - Si{kn)] 

+ cos{kr,) [Ci{[l + c]krs) - Ci{kr,)] - -^^^^H^ I , (9) 



(1 + c)krs 



^ http://camb.info/ 
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where Si and Ci are the sine and cosine integrals, respec- 
tively: 



Si{x) = 



smit) 
t 



dt, Ci{x) = — 



cos(f) 
t 



dt. 



(10) 



Following Bullock et al. (2001), we approximate the depen- 
dence of the concentration c on M and z as 

where M*(2:) is the characteristic mass scale at which 
i^{M,z) = 1; M,(z = 0) ~ 5 X 10^^ h'^ Mq. 

In the 1-halo term [eq. ((2])] we set s = 2, in analogy with 
the corresponding term for the dark matter power spectrum, 
if {Ngai{Nga,i — 1)) > 1. Otherwise we set s = 1 since if the 
halo contains only one galaxy, it will sit at the center. Taking 
into account that (A/'gai (iVgai - 1)) ^ 2{Nccn)(N,^t) + {N,^tf 
and that only the galaxies that are not at the center get 
factors of Ugai(fc, M) ~ iidm(fc, M) we have; 
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and pm/M — 3/(47r_R^). On large scales, where the 2- 
halo term dominates, itdm(fc,A/) — > 1 and P^^{k,z) ~ 
fcgai-Piin(fc, 2) with: 



fega.(.)= / duf{u)(P^)h{M,z)^-^ 

V M J ngai 



We also define the effective large-scale bias, bcs{z), as 
dn (Afgai) 



feoff (z) 



dM- 



dM rigai 



-&(M,2), 



and the effective mass of the halo, Mcb, 

MM^)= I dM^M^. 

Jm "g-i 



(15) 
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4 ANGULAR POWER SPECTRUM OF 
INTENSITY FLUCTUATIONS 

The angular power spectrum P{kg) of intensity fluctuations 
due to clustering of sources fainter than some flux density 
limit Slim is a projection of the spatial power spectrum of 
such sources at different redshifts z, Pgai(fe, z). In the Limber 
approximation, valid if the angular scale is not too large (i.e. 
27rfce ^ 10), the relation between Pgai(fc,z) and P{kg) is: 



Pike 



dz Pgi 



where dS/dz is the redshift distribution of the cumulative 
flux of sources with S ^ Sum 



^ = j''"^ d\og,,{S)SmS, z), z] ^, (19) 

<j)(L, z) is the epoch-dependent comoving luminosity func- 
tion per unit interval of logj^g(L), and dVc is the comoving 
volume element, dVc = X^dx, x{z) being the comoving ra- 
dial distance: 



X(^) = 



dz' 



Ho Jo y'nmoil + Z')3 + {1 ~ n^o) 



(20) 



Poisson fluctuations add a white noise contribution to the 
power spectrum of fluctuations: 

f-Sii„ 



dN 



rflogio S 



S dlogioS, 



with 
dN 



dz (j)[L{S, z), z 



dVc 

dz 
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We have computed the functions dS/ dz for each galaxy pop- 
ulation using the the cosmological model specified in §[T] 
and the evolutionary models briefiy described in §(2] As 
mentioned in §[2l the functions dS/dz for proto-spheroidal 
galaxies at frequencies ^ 353 GHz (A > 850 nm) have been 
scaled down by constant factors to comply with the mea- 
surements or the best estimates of the shot-noise levels. At 
higher frequencies our model accurately fits the observed 
source counts and therefore provides directly reliable esti- 
mates of the shot noise level. 

We have chosen to deal with the shot noise and the 
clustering contributions to the power spectrum of CIB fluc- 
tuations independently of each other because the former 
are independent of the parameters describing the cluster- 
ing and are strongly constrained by the available source 
counts. Moreover, when only relatively low resolution data 
are available, as is the case for Planck, there is a degeneracy 
between the shot-noise and the 1-halo clustering term. As 
clearly highlighted by Planck Collaboration (2011), an un- 
supervised least-square flt of the full CIB power spectrum 
measured by Planck, taking the shot-noise amplitude as a 
free parameter, leads to fits of similar quality with and with- 
out a substantial contribution from the 1-halo term. But fits 
with a low contribution from the 1-halo term imply shot 
noise amplitudes far in excess of those estimated from the 
source counts. The higher resolution of Herschel, SPT and 
ACT data breaks the degeneracy at !/ J5 600 GIfz and at 
y ^ 220 GHz, respectively, allowing a direct estimate of the 
shot-noise amplitude. 

As for the halo model, we have considered two dis- 
tinct populations, i.e. proto-spheroidal galaxies and late- 
type galaxies, both quiescent and starbursting. Taking into 
account the constraints on clustering of late-type galaxies 
coming from IRAS data (Mann et al. 1995; Hawkins et al. 
2001) we find that the contribution of these sources is al- 
ways sub-dominant and, correspondingly, their halo model 
parameters are very poorly constrained. Moreover the val- 
ues of Msat and cr(logjQ M) are poorly constrained also for 
proto-spheroidal galaxies (Planck Collaboration 2011). We 
have therefore fixed Msat ~ 20Mmin and cr(logio M) — 0.6 
[within the ranges found by Tinker & Wetzel (2010) from 
clustering studies of optical galaxies] for both populations, 
and Mmin.iato-tvDc = lo" Mp, and Qfsat.iatc-tvDc = 1. We are 
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Figure 1. CIB angular power spectra P{kg) at sub-mm wavelengths. Data from Planck Collaboration (2011) and i7ersc/iei/HerMES 
(Amblard et al. 2011). For the i/ersc/ieZ/HerMES data at 350 and 500 lira we have adopted the values corrected by Planck Collaboration 
(2011). At 250 we have used the values given by Amblard et al. (2011), that may be underestimated because of an over-subtraction 
of the cirrus contamination and a slight overestimate of the effective beam area. The conversion from the multipole number I used by 
Planck Collaboration (2011) and the wavenumber fc (arcmin"^) is fc = <?/(2 X 180 X 60). The lines show the contributions of the 1-halo 
and 2-halo terms for the two populations considered here [spiral and starburst (SS), and proto-spheroidal (PS) galaxies]. The magenta 
horizonal lines denote the shot noise level. 
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Figure 2. CIB angular power spectra P{kg) at mm wavelengths and 150 X 220 GHz cross spectrum. Data from Planck Collaboration 
(2011) and SPT (Hall et al. 2010; Shirokoff et al. 2011). The ACT data (Dunkley et al. 2011; Das et al. 2011) are in good agreement 
with the SPT ones and are not plotted to avoid over-crowding the figure. The lines have the same meaning as in Fig. [T] 
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then left with only 2 free parameters, i.e. Mmin and Osat, for 
proto-spheroidal galaxies. 

The angular power spectra of CIB anisotropies at 217, 
353, 545, and 857 GIfz on the multipole range 200 ^ ^ ^ 
2000 have been determined by Planck Collaboration (2011) 
using Planck maps of six regions of low Galactic dust emis- 
sion with a total area of 140 deg^ . In the same paper, the 
power spectrum measurements by Amblard et al. (2011), 
using i/ersc/ie//SPIRE data at 250, 350, and 500 /^m and 
extending down to sub-arcmin angular scales, i.e. up to 
£ 2 X 10*, were re-analyzed. It was found that Amblard 
et al. (2011) overestimated the correction for contamination 
by Galactic cirrus. Moreover, the diffuse-emission calibra- 
tion of SPIRE data was improved using the more accurate 
Planck/TIFI calibration. We have used the Amblard et al. 
(2011) data as corrected by Planck Collaboration (2011) at 
350 and 500 nm. No correction could be applied at 250 /im 
so that the data points at this wavelength could be under- 
estimated. 

Power spectrum measurements at mm wavelengths 
(around 150 and 220 GHz) have been obtained with the 
SPT and the ACT (HaU et al. 2010; Dunkley et al. 2011; 
Shirokoff et al. 2011; Das et al. 2011). The subtraction of 
the other components (CMB, Sunyaev-Zeldovich effect, ra- 
dio sources) has been done using the best fit values given 
in the papers. Note that the units quoted as fiK^ are ac- 
tually /iK^xsr. The conversion factor from these units to 
JyVsr is ~ [24:.8(x'^ / smh{x/2))f /[2£{£ + 1)]. The factor 



[24.8(a:7sinh(cc/2))] = 



1.55 X lO'^ at 150 GHz and 



2.34 X lO'^ at 217 GHz. 
The angular correlation function w{9) for a single source 
population writes, in terms of the 2D power spectrum P{ke): 



w(6) = 2-K I keP{ko)Jo{2Tvkee)dk0, 



(23) 



where Jo is the Bessel function of order 0. 

Here we have two sub-populations, proto-spheroidal and 
spirals+starburst galaxies, with different clustering proper- 
ties. If their cross-correlations can be ignored the signal for 
the whole is given by (Wilman et al. 2003): 



wtot{0) = fpswps{e) + fsswss{0), 



(24) 



where fps and fss are the fractional contributions of proto- 
spheroidal and spirals+starburst galaxies, respectively, to 
the total counts: 



PS/SS 



/ dzA/"ps/ss(z) 
/ dzJ\ftot{z) 



(25) 



Af being the redshift distribution. Ignoring the cross- 
correlations between the two source populations is justi- 
fied because of the widely different redshift distributions 
implied by the adopted evolutionary model: as mentioned 
in §[2l proto-spheroidal galaxies are associated to galactic- 
size halos virialized at z^ir ^ 1.5 while disk (and irregu- 

lar/starburst) galaxies are associated primarily to halos viri- 
alizing at z^ir 1.5. 

The spatial correlation function ^(r, z) is the Fourier 
anti-transform of the 3D power spectrum: 



1 

2^ 



sin(fcr) 
kr 



3x10' 



250 urn 
350 urn 
500 urn 



Proto-spheroids 




dk. 



(26) 



Redshift z 

Figure 3. Redshift evolution of the galaxy intensities at 250, 350, 
and 500 fira yielded by the model compared with the observation- 
based estimates by Amblard et al. (2011). The different line styles 
correspond to the different sub-populations, as specified in the 
inset. Some caution is needed in interpreting these data, in view 
of the problems pointed out by Planck Collaboration (2011; see 
text). 



The clustering radius ro{z) is defined by ^(ro, 



5 RESULTS 

Figures [1] and [2] compare the best fit model power spec- 
trum with Planck, Herschel, and SPT data in the wavelength 
range 250 /^m— 2 mm. The agreement is generally good ex- 
cept at 250 /im where the model is consistently above the 
data points by Amblard et al. (2011) which, however, could 
be underestimated (see §|4]). As mentioned above, we have 
only two free parameters, i.e. the minimum mass and the 
power-law index of the mean occupation function of satellite 
galaxies of proto-spheroidal galaxies. The constraints we ob- 
tain are log(Afmin/MQ) = 12.24±0.06 and Qsat = 1.81±0.04 
(Icr). The nominal errors on each parameter have been 
computed marginalizing on the other and correspond to 
Ax^ = 1. We caution that the true uncertainties are likely 
substantially higher than the nominal values, both because 
the model relies on simplifying assumptions that may make 
it too rigid and because of possible systematics affecting the 
data. 

The fact that the same values of these parameters ac- 
count for the clustering data from 2 mm to 250 /im con- 
firms the conclusion by Planck Collaboration (2011) that 
CIB fluctuations over this wavelength range are dominated 
by a single sub-population of dusty galaxies. According to 
our model, this sub-population is made of proto-spheroidal 
galaxies making most of their stars at z > 1. We find that 
only at 250 /im other dusty galaxy populations, normal disk 
and starburst galaxies, make a significant, but still sub- 
dominant contribution to the clustering power spectrum. As 
shown by Lapi et al. (2011), according to our model, proto- 
spheroidal galaxies also account for the bulk of the CIB in- 
tensity in this wavelength range, consistent with the finding 
by Planck Collaboration (2011) that the CIB anisotropies 
have the same frequency spectrum as the CIB intensity. 
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Figure 4. Angular correlation function of sub-mm galaxies. Data from Cooray et al. (2010) and Maddox et al. (2010). In the left-hand 
panel the dashed and dash-dotted lines show the contributions of PS and SS populations, while the solid line denotes the total angular 
correlation function. In other two panels, the dashed and solid lines show the global model correlation functions for the limiting flux 
densities adopted in the analyses by Maddox et al. and Cooray et al., respectively. The horizontal dotted line shows the zero level. 



Our estimate of the minimum mass is higher than, but 
consistent, within the errors, with those found by Amblard 
et al. (2011) considering a single galaxy population and 
5 free parameters per frequency (but one of the parame- 
ters is unconstrained by the data within the prior range): 
log(Mn,in/M0) = 11.1+0^ at 250 /im, log(M^in/M0) = 
11.51^:^ at 350 Atm, and log(Mmi„/Mo) = H-Slo.s at 
500 /im. Our value of asat is also consistent with those by 
Amblard et al.: Osat ~ 1-6!'^q'2 at 250 /.im, Qsat = I.S^q'^ 
at 350 /im and 500 /im. In Planck Collaboration (2011) two 
or three free parameters per frequency were used; the de- 
rived minimum masses are in the range log(Minin/Af0) = 
11.8 - 12.5. 

There is however an interesting difference with Planck 
Collaboration (2011), due to the different redshift distri- 
butions of sources. The crossover between the 1-halo and 
the 2-halo term occurs, according to the model by Planck 
Collaboration (2011), at multipole numbers ranging from 
£ ~ 800 at 857 GHz (350 /tm) to ^ ~ 1200 at 217 GHz 
(1.38 mm) corresponding to angular scales ranging from 
e ~ 180 X 60/£ ~ 13.5' at 857 GHz to 9' at 217 GHz. Accord- 
ing to the Bethermin et al. (2011) model used in that paper, 
the contribution to the CIB intensity at 857 GHz peaks at 
z ~ 1 where the angular scale of 13.5' corresponds to a phys- 
ical linear scale L ~ 6.5 Mpc; at 217 GHz the bulk of the CIB 
contribution comes from z > 2 where an angular scale of 9' 
corresponds to a physical linear scale L ~ 4.5 Mpc. The non- 
linear masses corresponding to an overdensity Ac = 1.68 on 
these scales are Mni{z = l;L = 6.5Mpc) ~ 7 x lO" Mq and 
Mni(z = 2;L = 4.5Mpc) ~ 8 x IO^^A/q, respectively. For 
comparison, the characteristic non-linear masses computed 
from <j{Mt,,z) = 1.68 ((T(Af, ,2) being the rms overdensity) 
areM.(2 = 1) = 2xl0" Mq and M*(2 = 2) = 7.3x10® Mq. 
This suggests that structures going non-linear on the con- 
sidered scales are extremely rare at the corresponding red- 
shifts. This potential difficulty is eased in our model because 
the crossover scales are lower by almost a factor of 2. As 
shown by Figs. [T] and O the l-halo/2-halo crossover occurs 
at ^ ~ 1450 at 857 GHz and at ^ ~ 2100 at 217 GHz, corre- 
sponding to angular scales of 7.4' and 5.1', respectively. 

Our value of i\fmin implies an effective halo mass 
[eq. p7p ] at z ~ 2 of proto-spheroidal galaxies, making up 
most of the CIB, M^fi ~ 5 x lO^'^ Mq, close to the estimated 
halo mass of the most effective star formers in the universe. 



Tacconi et al. (2008) estimated their mean comoving density 
at 2 ~ 2 to be ~ 2 X 10"* Mpc"^. For the standard ACDM 
cosmology this implies that they are hosted by dark matter 
halos of ~ 3.5 x 10^^ Mq (Dekel et al. 2009). 

Figure [3] compares the flux density coming from differ- 
ent redshifts, dS/dz [eq. (|19|) ]. predicted by the model at the 
SPIRE wavelengths with the best fit estimates by Amblard 
et al. (2011). Planck Collaboration (2011) give (their Ta- 
ble 7) the best fit values of the redshift-independent volume 
emissivity, jeff, for z > 3.5 [their eq. (43)]. The values of jcs 
given by our model (52, 175, 265, and 205 Jy sr~^ Mpc^^ at 
217, 353, 545, and 857 GHz, respectively) are consistent with 
the best-fit results. 

As for the angular correlation function, 1^(9), of de- 
tected SPIRE galaxies, Cooray et al. (2010) reported mea- 
surements of w{6) for sources brighter than 30 mjy at all 
SPIRE wavelengths and inferred values of log(Afmin/MQ) 
ranging from 12.6^0.6 at 250 /im to 13.5t ".o at 500 /im. On 
the other hand, Maddox et al. (2010) did not detect a sig- 
nificant clustering for their 250 /im selected sample with a 
flux limit of 33mJybeam~^, but detected strong cluster- 
ing at 350 /im and 500 /im, albeit with relatively large un- 
certainties. Our model entails a relationship between the 
far-IR/sub-mm luminosity of proto-spheroidal galaxies (that 
provide the dominant contribution to w{6), see the left-hand 
panel of Fig. |4]) and the associated halo masses (Lapi et al. 
2011). For the flux density limit adopted by Cooray et al. 
(2010), 30 mJy at all SPIRE wavelengths, the model yields 
log(Mmin/MQ) ~ 12.3, 13, and 13.2 at 250, 350, and 500 /im, 
respectively, while for the flux density limits of Maddox et al. 
(2010; 33, 36, and 45 mJy) we have log(Mmin/MQ) ~ 12.3, 
13.1, and 13.4. The corresponding predictions for w{6) are 
compared with the data in Fig. [l] The agreement of the 
model with the data is generally good, although the situ- 
ation at 250 /im is unclear since there is a discrepancy be- 
tween the Cooray et al. (2010) and the Maddox et al. (2010) 
results. 



6 DISCUSSION AND CONCLUSIONS 

According to the Granato et al. (2004) model, the steep 
portion of sub-mm counts is dominated by massive proto- 
spheroidal galaxies in the process of forming most of their 
stars on a timescale varying with halo mass (shorter for more 
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massive galaxies), but typically of ~ 0.7 Gyr, i.e. with a duty 
cycle of ~ 0.2 at z ~ 2, where their redshift distribution 
peaks. As shown above, this model allows us to reproduce 
the power spectrum of GIB fluctuations over a broad fre- 
quency range, from 250 /im to a few mm, with only 2 free 
parameters. The model also yields an effective volume emis- 
sivity at different redshifts consistent with observational es- 
timates. The derived effective halo mass, Moff ~ 5x10^^ Mq, 
is close to that estimated for the most efficient star-formers 
at 2 ~ 2. 

The multipole number at which the 1-halo term starts 
exceeding the 2-halo contribution to the clustering power 
spectrum ranges from i — 1450 at 857 GHz to £ = 2f00 at 
217 GHz. These values are almost a factor of 2 higher (and, 
consequently, the corresponding angular scales are almost 
a factor of 2 lower) than those found by Planck GoUabora- 
tion (2011). Since, at the redshifts where the contribution 
to the GIB intensity peaks, the corresponding masses are 
well above M« , this difference translates into a much larger 
abundance of the relevant halos. 

Alternative models make quite different predictions for 
the clustering properties of sub-mm galaxies. A widespread 
view is that these objects are powered by major merger 
events. Two major theories have been worked out in this gen- 
eral framework. One view is that sub-mm galaxies are mas- 
sive objects, seen during a short-duration, intense, merger- 
induced burst of star formation (e.g. Narayanan et al. 2009). 
Since massive galaxies are rare at high-z and, because of the 
short duration of the burst, only a small fraction of them are 
in the sub-mm bright phase at a given time, this scenario has 
difficulty in reproducing the observed counts. This difficulty 
may be overcome assuming an extremely top-heavy initial 
stellar mass function that would allow much less massive 
(hence far more abundant) galaxies to reach the required 
luminosities (e.g. Baugh et al. 2005; Lacey et al. 2010). 

The clustering implied by the latter scenario has been 
investigated by Almeida et al. (2011) who found, at z = 2, 
a comoving correlation length of ro — 5.6 ± 0.9/i~^Mpc 
for galaxies with 850 fim flux densities brighter than 5 mjy 
or an effective bias factor &cft = 2.3; for galaxies with 
5*450^111 > 5mJy they found b^s = 2.1. Our model implies 
log(Mmin/Af0) ~ 12.4 for sources with S450;jm > 5mJy and 
log(Mniin/M0) ~ 13.0 for sources with S'ssoMm > 5mJy. The 
corresponding values of the clustering radius and of the effec- 
tive bias factor are ro ~ 11.2 Mpc, bcs = 4.3 at 850 /im, 
and ro — 7.3 Mpc, &eff = 3.1 at 450 /xm. The study by 
Kim et al. (2011) confirms that the clustering data require a 
higher amplitude of the 2-halo term, i.e. more massive halos 
than implied by the major mergers plus top-heavy initial 
stellar mass function scenario. 

Dave et al. (2010) investigated the clustering properties 
of rapidly star-forming galaxies at z ~ 2 in the framework 
of a very different scenario based on cosmological hydrody- 
namic simulations whereby the star formation is not pow- 
ered by mergers but by steady gas accretion and cooling 
that can fuel the star formation for several Gyrs. In this 
scenario typical sub-mm galaxies at z = 2 live in massive 
(~ 10^^ Mq) halos and have a duty cycle ~ 50%. They are 
expected to be strongly clustered, with a clustering radius 
ro ~ 10/i~^Mpc and a bias factor of ~ 6. These values 
are well in excess of those following from our analysis which 



yields, for the bulk of galaxies at z ~ 2, ro — 6.9 h ^ Mpc, 
&cff ^ 3. 

These results illustrate the power of accurate measure- 
ments of the GIB power spectrum and of the correlation 
function of galajcies at (sub-)millimeter wavelengths to dis- 
criminate among competing evolutionary models for the 
population of dusty galaxies. 
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